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Abstract 



We propose a new approach to construct QCD sum rules for the irNN cou- 

oo ■ 

■ pling constant, g, starting from the vacuum-to-pion correlation function of 

o : 

the interpolating fields of two nucleons and taking its matrix element with 
respect to nucleon spinors. The new approach with the projected correlation 
function is advantageous because even in the chiral limit the dispersion in- 
tegral can be parametrized with well-defined physical parameters. Another 

x 

5_i ■ advantage of the new approach is that unwanted pole contribution is projected 

out. Calculating the Wilson coefficients of the operator product expansion of 
the correlation function up to O \ Mb~j and O (mj where Mb and are 
the Borel mass and the pion mass, respectively, we construct new QCD sum 
rules for the irNN coupling constant from the projected correlation function 
with consistently including 0(m n ) corrections. By numerically analyzing the 
obtained four sum rules we identified the most prominent one. After roughly 
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estimating errors we obtaind, g = 10 ± 3, as a result of this sum rule, which 

is in reasonable agreement with the empirical value. It is also found that the 

0(m n ) correction is about 5%. 

PACS number(s): 13.75.Gx, 11.55.Hx, 24.85.+p 
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I. INTRODUCTION 



Within the framework of the QCD sum rule proposed by Shifman, Vainshtein and Za- 
kharov ffl, much work has been done for the meson-baryon-baryon coupling constants 0-0. 
Sum rules have been constructed in two formulations: one is to start from the vacuum-to- 
vacuum matrix element of the correlation function of the interpolating fields of two baryons 
and one meson and the other is to start from the vacuum-to-pion matrix element of the 
correlation function of the interpolating fields of two baryons. In the pionnering work by 
Reinders, Rubinstein and Yazaki the former approach was taken but in the following works 
including the one by the same authors the latter approach was employed since some diffi- 
culties seem to exist in the former approach |§. 

Therefore, we discuss only the formulation with the vacuum-to-pion matrix element of 
the correlation function of the interpolating fields of two nucleons: 

U{p,k) = -i f d 4 xe tpx (0\T{ri(x)rj(0)}\ir(k)), (1) 

where rj(x) is the interpolating field for the nucleon, |vr(fc)) is the pion state with momentum 
k normalized as (ir(k')\7r(k)) = 2fco(27r) 3 5 3 (fc' — k) and the isospin is neglected for simplicity. 
The correlation function can be split into four Dirac structures, 275, ij^fi, ijsfk and 75a ' jjLV p tl k v ' . 
The coefficients of the Dirac structures are functions of Lorentz invariants, p 2 and pk, and are 
called invariant correlation functions. Reinders, Rubinstein and Yazaki constructed a sum 
rule for the pion-nucleon-nucleon (jrNN) coupling constant in the chiral limit with vanishing 
pion momentum from the invariant correlation function with the Dirac structure 275 in the 
leading order of the operator product expansion (OPE). They showed that if the sum rule is 
divided by the odd-dimensional sum rule for the nucleon mass, the result is consistent with 
the Goldberger-Treiman relation with gA — 1, where qa is the axial charge of the nucleon ||. 
Shiomi and Hatsuda improved the sum rule by taking into account higher dimensional terms 
and a s corrections of the OPE and also the continuum contributions {3|. They showed that 
even after these corrections are taken into account, the Goldberger-Treiman relation with 
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qa = 1 holds as long as the same continuum threshold is taken in the sum rule for the wNN 
coupling constant and the odd-dimensional sum rule for the nucleon mass. In Ref. || Birse 
and Krippa pointed out that in the chiral limit the vacuum-to-pion correlation function is 
obtained just by chirally rotating the vacuum-to- vacuum correlation function and therefore 
it is obvious that the ratio of the two sum rules is consistent with the Goldberger-Treiman 
relation, which is just a consequence of the chiral symmetry. Then, they tried to obtain a 
new sum rule for the irNN coupling constant which is not just a consequence of the chiral 
symmetry. They considered the invariant correlation function with the structure i^f^f, and 
obtained a sum rule by taking the limit, k — 0, after removing ijsfk. In all of these works 
the sum rule is constructed by making an "ansatz "for the absorptive part of the correlation 
function based on the effective lagrangian with the pseudoscalar coupling scheme. In Ref. 
Kim, Lee and Oka examined how the choices of the effective lagrangian, i.e. the pseudoscalar 
and pseudovector coupling schemes make differences in the sum rules for all the four Dirac 
structures. They concluded that only the invariant correlation function with the Dirac 
structure Jscr^p^k" , is independent of the two coupling schemes and obtained a sum rule 
by taking the limit, k — 0, after removing ^scr ^p 11 k u . They also claimed that their sum rule 
suffers from less uncertainties due to QCD parameters. 

In a recent work, we have studied in detail the physical content of the invariant corre- 
lation functions without referring to the effective lagrangian ||. We have shown that the 
coefficients of the double poles are proportional to the nNN coupling constant but that the 
coefficients of the single poles are not determined by the nNN coupling constant. In the 
chiral limit the double-pole terms survive only for the Dirac structures i^sfc and ^^c^p^ky . 
For these structures the single-pole terms as well as the continuum terms turn out to be 
ill defined in the dispersion integral. Therefore, the use of naive QCD sum rules for these 
structures in the chiral limit is not justified. 

In this paper we propose a new approach, a projrcted correlation function approach. We 
investigate the QCD sum rule for the nNN coupling constant derived from a projrcted cor- 
relation function by applying the procedure of Ref. [|T0 1 . An advantage of the new approach 



is that the discontinuity of the projected correlation function does not include ill-defined 
terms even in the chiral limit in contrast to the previous approaches. Another point is that 
we obtain sum rules not only in the chiral limit but also with consistently including 0{m- K ) 
corrections. 

This paper is organized as follows. In Sec. II, we summarize the physical content of the 
correlation function without referring to the effective lagrangian. We introduce the projected 
correlation function and study its physical content in Sec. III. In Sec. IV we construct QCD 
sum rules for the projected correlation function and investigate the numerical results of the 
obtained sum rules. Section V is devoted to summary. 



II. PHYSICAL CONTENT OF THE CORRELATION FUNCTION 

In this section, we summarize the physical content of the correlation function defined 
by Eq. (H) without referring to the effective theory ||]. The correlation function IT(j>, k) in 
Eq. ([I]) has a pole at p 2 = M 2 and (p — k) 2 = M 2 where p and p — k, respectively, become 
on-shell momenta for the nucleon where M denotes the nucleon mass. The ttNN coupling 
constant, g, is defined through the coefficient of the pole as 

u{pr)(j> - M)U(p, k)($-fc- M)u(qs)\ p 2 =M 2 i(p _ k) 2 =M 2 = -i\ 2 gu{pr)j 5 u(qs) , (2) 

where q = p — k, u(pr) is a Dirac spinor with momentum p, spin r and is normalized as 
u(pr)u(pr) = 2M. A is the coupling strength of the interpolating field to the nucleon. 

In addition to the pole singularity, H(p, k) has a branch cut singularity starting from the 
threshold of the pion-nucleon channel, p 2 = (M + m^) 2 , to infinity where m n is the pion 
mass. We restrict ourselves to the positive energy region, p > 0, for simplicity. 

In order to classify the singularity of U(p,k), we define the ttNN vertex function, 
T(p,q',k f ) (p — q' + k'), and the pion-nucleon T-matrix, T(q' , k! , q, k) (</ + k! = q + k), 
by 



T(p,q',k') = ^-M) 



i / d 4 xe ipx (0\T[ri{x)f]{0)]\ir{k')) 



f-M), 



T{q\ k', g, k) = d' - M) 



The correlation function is related to the vertex function as 



i rf 4 xe Vx (7r(A; / )|T^(x)r7(0)]|7r(A;)) (j-M). (3) 



p 2-M 2 \p-k) 2 -M 2 ' 

Let us regard H(p, k) as a function of the center-of-mass energy, i.e. po in the frame p = 0, 
and consider the discontinuity of the correlation function. We adopt the following notation 
for the dispersive (continuous) part and the absorptive (discontinuous) part, respectively: 

ReF(p) = - [F(p)\ p o =p0+iT1 + F(p)\ p0=p o_ irj ] , 

ImF(p) = ~r [F(p) \ p o =p0+iv - F(p)\ p0=p o„ iri } . (4) 
The absorptive part of the correlation function, Imll(p, k), can be written as 
Iml% k) = H,+ M){lm (p2 _ M2)( ^_ k)2 _ M2) ReT(p, q, k) 

_ M 2 )((p — kf ~ M*) ""^' 9 ' k) }<> ~ * + M) - (5) 

The first term represents the pole contribution and the second term represents the continuum 
contribution. 

When the center-of-mass energy, p , is above the threshold of the pion-nuclon channel, 
M + m n , but below that of the next channel, only the pion-nucleon channel contributes in 
the intermediate states and the absorptive part of the vertex function, Imr(p, q, k), is given 
by 

1mT(p,q,k) = -~ I ^T(p,q',k')(4' + M)T*(q',k',q,k), (6) 



where t is defined through po = V M 2 + t 2 + ym 2 +t 2 . q' and k! are the on-shell momenta 
for the nucleon and the pion, respectively. 

From Eq. (|(J), one sees that in the vicinity of the pion-nucleon threshold T can be ex- 
panded as 
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Imr(p, q, k) 



ReT(p,q,k) 



(po < M + m*) 

Git + 0(t 3 ) (po >M + m 7r ) ' 

(7) 

G + G x t + 0(r 3 ) (p <M + m 7r ) 



G + O(t 2 ) (p >M + m 7r ) 

where r = it. In Eq. ([7]) the behaviour of the dispersive part is determined from that of the 
absorptive part by analytic continuation. From Eq. ([I]) one sees that 

oo (m n 7^ 0) 

P o=M+m 7r -e I finite constant (m„ = 0) 



d 

- — Rer(p, q, k) 



d 

—ReT(p,q,k) 



= 0. (8) 

p =M+m w +e 



We found in Ref. that in the chiral limit the double-pole terms survive only for the 
invariant correlation function with the Dirac structures 275^; and 75a ^pf 1 k v . The coefficients 
of the single poles for the structures 275^ and 75a i]LV p p 'k v involve the derivatives of the vertex 
functions with respect to p 2 and are therefore ill defined due to Eq. (f|). Therefore, the use 
of naive QCD sum rules for these structures in the chiral limit is not justified. 

A reasonable procedure to make an ansatz for the absorptive part of the invariant cor- 
relation function to be used in the sum rule would be first to transform it, for instance by 
multiplying by p 2 — M 2 , to a form which does not include ill-defined terms. However as 
pionted out in Refs.||U one is unable to find a reasonable result using this proceedure. In 
this paper, we will propose a different way to construct a sum rule. 



III. PHYSICAL CONTENT OF THE PROJECTED CORRELATION FUNCTION 

We define the projected correlation function, tcNN vertex function and pion-nucleon 
T-matrix by 

n+Qor, qs, k) = u(pr)j ^-{p, k)'jou(qs), 
T + (pr, qs, k) = u(pr)T(p, q, k)u(qs), 

T + (qs, k, q's', k') = u(qs)T(q, k, q', k')u(q's'). (9) 
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The projected vertex function is related to the projected correlation function as 



T + (pr, qs, k) = (p - M)(p - E q - u k )Ii+(pr, qs, k), 



(10) 

where E q = y/M' 2 + q 2 and u k = \Jm 2 + k 2 . It should be noted that H + (pr, qs, k) has poles 
at po = M and po = E q + u k but not at po = —M and po = —E q + u k . Then the absorptive 
part of the projected correlation function, Imil + , can be written as 

t i r / / \ xf A/f \ Rer + (pr, qs, k) ReT + (pr,qs,k) 

han+ipr, qs, k) = 7i5{p - M)-— — - 7r5(p -E q - w fc -— — 

E q + uo k — M E q + uo k — M 

lmT + (pr,qs,k). (11) 



+Re- 



(po - M) (p -E q -u k 

Let us now consider the last term of Eq. (|TT|). When the center-of-mass energy, po, is above 
the threshold of the pion-nucleon channel, M + m^, but below that of the next channel, 
Imr + (pr, qs, k) is given by 

lmT + (pr,qs, k) = -7—^ / 77-^ r +(P r > l' s ^ k')T*(q's', k', qs, k) 

= -j-^ 2 g(Po,t 2 ) I ^u{pr)i lh {i' + M)T*(q>k',q,k)u(qs), (12) 

where q' and kl are on-shell momenta for the nucleon and the pion, respectively, and g(po, t 2 ) 
is defined by 



T+(pr, q s , k') = X u(pr)i<y 5 u(q's')g(po,t ). 



(13) 



Since w(p?")^75 W + M) — > u(pr)i r y 5 (—t 2 /2M — 7^) — > as p — > M + m, Y + behaves in the 
vicinity of the pion-nucleon threshold as 



Imr + (pr, qs, k) 



ReT + (pr, qs, k) 



(po < M + nv) 

0(t 3 ) (p >M + m 7r ) 
G + + 0(r 3 ) (p <M + m, 
G + + 0(t 2 ) (p >M + m, 



(14) 



which is in contrast to Eq. (|7]). 

Here we should discuss the behaviour of Imn + in the chiral limit with k = 0: We obtain 



S 



d 

ImII + (]?r, qs, k) = —tt5'(po — M)ReT + (pr, qs, k) + n5(po — M)— — ReT + (pr, qs, k) 

opo 

+Re (p - M )2 Imr +(P r 'ff g ' fc )- ( 15 ) 

The coefficient of the single pole term vanishes since ^Rer + behave as (po— M) 2 for po < M 
and po — M for po > M. The last term behaves as 8(po — M)(po — M) since Imr + behaves 
as 6(po — M)(p — M) 3 . Therefore, the second term and third term on the right-hand side 
of Eq. flTSp are well defined in the dispersion integral. 



IV. CONSTRUCTION OF SUM RULES 

For the sake of completeness we first summarize the derivation of the Borel sum rule 
for the projected correlation function. The dispersion relation for the projected correlation 
function in the variable po is 

n+ p ,P = / dp ——. 16 

J Pa-Po+iV 

By splitting the projected correlation function into the even part, n +even = |[n(po,p) + 
n(-p ,p)], and odd one, Il +odd = ^[U(p ,p) -U(-p ,p)), we rewrite Eq. (|l|) in terms of 
the even and odd parts as 



n+even (p%,p) = ~~ f d Po 2^ ,2 ImII + (Po» P) ' 
71 J pi- ptf 

If 1 

n +0 dd(Po>P) = / d Po— ^ Imn +(Po,p). (17) 

n J pi — Pn 



Applying the Borel transformation, 

(pir ( d\ n 



j b = lim , , 

-Ipoo (n-1)! V dpi 

-p2/n=M B 2 



where Mb is called the Borel mass, we get 

If p' ( p' 2 \ 
L B [U +evea (pl,p)] = - J dp' j^exp [-jf^j Imn + (p[,,p), 

■ 2 



If 1 ( P \ 

L B [U +odd (p 2 ,p)] = - J dp' j^exp f -j^J Imn + (p(,,p). (18) 



These equations with the correlation functions on the left-hand side evaluated by the OPE 
are the Borel sum rules. 

Now, we calculate the OPE of the two-point correlation function. We consider the 
following correlation function with charged pion, 

U(p,k) = -iJ d'x^iOlT^Um^ik))- (19) 



with the interpolating field for the proton and the neutron as [111 



Vn = -e & c [daC , 7 / A]757 M w c , 

where u and d are the up-quark field and the down-quark field, respectively, C is the charge 
conjugation operator and a, b and c are color indices. The Wilson coefficients for eq.(H5) 



have been calculated in Ref. 0(k°) terms for the structure 275, and 75a fMU k fl p lJ 

are obtained in Ref. [^,[| , [[J and ||, respectively. Ref || also gives O(k) terms for the 
structure ^75. Here, we also need 0{k) terms for the structure 275^, 275^ and iiv k p "p v . We 
calculate the Wilson coefficients of the short- distance expansion in two steps: we perform 
the light-cone expansion of the correlation function first and the short- distance expansion 
of the light-cone operators second. The reason for doing this is to use the parameterization 



of the vacuum-to-pion matrix elements of the light-cone operators given in Ref. [12 



(0| d(0)i l5 u(x)\rr + (k)) = V2 fir ™" (l - % -kx + ((kx) 2 )) , 

m u + m d V 2 ' J 

(0\d(0) l5 a^u(x)\K + (k)) = V2i(k^ - k v x u ) (l - % -kx + O ((kx) 2 )) 

(0\d(0) l5llM u(x)\7i + (k)) = V2ikJ w (l - % -kx + ((kx) 



+ ^2 tkfiX ^^hL ^ - l -k x + ((kx) 2 
fJ 2 ' 



-V2i Xfl kx^— [l-~kx + ((kx) 2 ) 
18 V 2 ^ ' 

(0|d(0)5f s G Q/3 (0)7 M u(a;)|7r + (A;)) = -y/2i(k^g a ^ - k^^f^S 2 



1 .. 7 + 6e . , w ., , 



x ikx h O ( (kx 

V3 84 V v 

+V2{k a xp - kpx^k^fj 2 
10 



(20) 
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f n is the pion decay constant, m u (m^) is the up(down)-quark mass, G p(T is the gluon field 
tensor and G a p = ^Sa^G^ ' . 5 2 is defined through g s (0\d~G ^ v 1 u u\tx + (k)) = \[2%f. K b 2 k lL where 
g s is the coupling constant of the QCD. e is a parameter associated with the deviation of 



twist-4 pion wave functions from their asymptotic forms [13 



Going through the above two steps, we obtain the OPE of the correlation function as 



1% k) 



V2 



«75 



4tt 2 k 

+^75 



p 2 ln(— p 2 ) 



(qq) 

f* ' P 



1 



rr 2 (qq)(^G 2 Y 



6 



p 2 H-p 2 )2f n + H-p 2 )2fj 2 + 



-l^^k u 
1 



7T 



{m)9s{q(T^G^q)' 



ln(— p 



2 (qq) 1 16tt 2 



3/. + P 2 3 



P 



fn(qq) 



ln 2 (qq) 2 tt 2 a s 2 



+ - 



2tT 



Ugs(q^G^q) 



+i-y 5 pk 

+i>y 5 fipk 

+i'y 5 ]tpk 



ln(— p 2 



(qqY 



U J + P 4 \ 6 A , 
2 \ 1 16tt 2 (gg) 2 1 4vr 2 (gg) 



* 2 (m)(^G 2 ) 1047T- 
54 

' 7r 2 (gg>(fG 2 )' 



27 



+ 



g s {qOu V G^q) 



+l5cr f iuP^k ,/ pk 



P 4 9 A > 6 9 A 
ln(V)(-2 / .) + i(-^) + i^(^) 



3 A 



+ 



P 



p u 



U8 + 32 % 2 fJ 2 (qq)- n - { ^G 2 i 



(21) 



63 " Jn ~ 9/^71 
where (O) denotes the vacuum-to-vacuum matrix element of the operator O, (qq) = (uu) = 
(dd) and G 2 = 2tr(G^G^). 

Some comments are in order here. The matrix elements of the four-quark operators and 
the mixed quark-gluon operators are approximated to factorize. Eq. ( pT|) is expressed in 
terms of the quark condensate, (qq), by the use of the Gell-Mann-Oakes-Renner relation, 
f 2 m 2 = —{m u + rrid)(qq). Though one might think that the terms with ry 5 ^ and r y^,a liu p^k v 
are of higher order in k than those with ^75 and 275^, the former and the latter are of the 
same order when sandwiched by u(pr) , yo an d 7o"u(qrs) in the center-of-mass frame: 
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u(pr)j ij 5 fcj u(qs) = (E q + M — u k )u(pr)^i^ Q u(qs) 



u k -E q -M 



Po 



u(pr)7o-f 5 a flu p> 1 k u -f u(qs) = -p {E q + M)u(pr)j ij 5 j u(qs) 

= (E q + M)u(pr)j i'j 5 fry u(qs). 

Thus, in Eq. we showed the terms which survive up to O (po 4 ) an d O (io k ) after pro- 
jection, and therefore the dimension of the operator does not play the role of the expansion 
parameter. 

Next, we express the right-hand side of the Borel sum rule, Eq. fll8|) , by physical quan- 
tities. We would like to make an ansatz for the absorptive part of the projected correlation 
function of Eq. |J) to be used in the sum rule. Since Eq. ( fL5|) does not include irregular 
terms, we make the following ansatz for the absorptive part of the projected correlation 
function for the 7r + -neutron-proton vertex: 

t rr i M a< A/r sR e r+(pr,qs,k) ReT + (pr,qs,k) 

Imn+(pr, qs, k) = ir6(p - M — — — - tt5{ Po - E q - u k )—— — 

E q + uj k — M E q + uo k — M 

+ [9(p - uo n ) + 9(-po - Imn° PE (pr, qs, k) 
= 5(p - M)u(pr)t l5 u(qs)V27T\ 2 j( M > k > 

-8(p -E q - u k )u(pr^Mqs)V2n\ 2 9 l\ + ^ ^ 

E q + u k - M 

+ [6(p - co w ) + 9(- Po - cu w )} Imn° PE (pr, qs, k), (22) 

where uj w is the effective continuum threshold of 11N or ttN channel and v2 on the right-hand 
side is the isospin factor. 

As is explained in Eq. (0), the 7rNN coupling constant is defined at p 2 = M 2 , (p — k) 2 = 
M 2 , i.e., 

9 = 9{ P o = M,k 2 = -m 2 n + ^ ] j, (23) 



while g(M, k ) and g(E q + uj k , k ) in Eq. Q2"2] ) refer to different kinematical points. However, 
if one expands them in m n taking k = one finds: 
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g(M,0) = g + O(ml), 

g(M + m^O) = g(M,0) + m7r 



dg(po, 0) 



dp 



O ( m l) = 9 + g'm* + O (ml) . 



Po =M 



Therefore, if we are interested only in the result of 0(m n ), the difference of g(M, 0) and g can 
be ignored. Then, substituting Eq. ([22) into the right-han side of Eq. ( |T8"D and expanding 
in 0(m n ) we obtain 

■ 2 \ 



b 2 [ M B \ 
2M ( 



ImIT + (pV, qs, k) 



-u(pr)ij 5 u(qs)^^ \ </ 



9'M + g 
2M 2 



M B 2 



V2\ 2 



M 2 \ ,„ 

up +(Cont ' 



Po 



M B Z 



Imn + (p'r, qs, k) 



2M 



-u(pr)iy 5 u(qs)j^ \ g - ^jMb 



V2X 2 



M 2 



Mf 2 



(Cont.), 



(24) 



i b \ 2M 1 ivib \ J-VJ-B 

where 0(m 2 ) terms are ignored. In Eq. (|2~4| ) (Cont.) denoes the continuum contribution 
coming from the last term of Eq. (^), which can be taken into account by multiplying the 
term p 2 ( n_1 ) ln(— p 2 ) in the OPE by 



C n (u) 



after Borel transformation. 



E 



U) 



n— 1 



x(k — 1)! \M E 



or 



exp 



M, 



2 ' 



B 



Substituting Eq. ( plj) and Eq. (|24|) into the left-hand side and the right-hand side of 



Eq. (|18D, respectively, we obtain sum rules for the coupling constant, 

V2X 2 ( M 2 \ 



2M ( g'M + g 
2 \ 9-- ~ 



M B 
Ait 2 } 



2M 2 
M B 2 C 2 M 



Mb 2 I tt^ exp 



M, 



B 



M, 



B 



(gg) 



+ (E q + M- u k )2f n 



C^KEg + M-u^fJ 2 



M, 



B 



vr 2 (qq)(^G 2 ) 



6 



n 2 ( (qq)' 



+ (E q + M — u k 



-u k - 18-^- + fA^G 2 } + (E q + M)16U(qq) 



M B 4 



— (Eg + M — „ fc )^^M! _ U f 



9 U 



7T 



-u k (E q + M)'—UUml(qq) 



448 + 32e 
21 



fJ 2 (qq) 



1 (qq) ,01 

3 U vr 



S G 2 ) 



13 



(25) 



and 



2M / 



2M 3 



V2X 2 



M 2 



M B 2 



V2 
4tt 2 ' 



- (/•;, • M)M_ 



3/. 



'(gg) 



+ (£, + M)2/ 7r 



1 



M ? 4 



7T 2 



+ M) 



16tt 2 „ , , 16 



71 



-(E q + M)'-\2Uml{qq) 



1 (qq)(^G 2 



fmqq) 



i{qq){^G 2 ) 



3 V 2 U 
= ^ s [ll +odd J, 



+ (E q + M)f 7T (-G 2 ) 

71 



(26) 



where we have used the parametrization, g s (qcr tlv G t11 ' 'q) = m^qq). 

From the sum rule for the vacuum-to-vacuum matrix element of the same operator, the 
strength A 2 is given up to the dimension seven by 



M o exp 



M 2 



M. 2 



b I 4tt 2 



M B 2 C 2 (u )(-(qq)) + 



* 2 (qq)(fG 2 } 

6M R 2 



OPE-Cont. 
even 



(27) 



and 



M 2 



Mf> 2 



4tH 



M 4 
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OPE-Cont. 1 



Qs G 2 ) + ^ ^ 2 



3M R 2 



(28) 



where Uq is the effective continuum threshold. 

Now, we eliminate A 2 by dividing the vacuum-to-pion sum rule by the vacuum-to- vacuum 
sum rule. Since we have even and odd equations for the vacuum-to-pion and the vacuum- 
to-vacuum correlation functions, respectively, we obtain the following four sum rules 
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Cont. 



9 2M2 Mb " 4M2 L B [n °P d E d - Cont -]' 

n 9' M2 _ V2M B " L B [UZ E d - d C ° nt -] (2q) 
9 2M " AM L B [U^J d - Cont -y 1 ' 

The procedure of extracting the coupling constant goes as follows. We regard the right- 
hand side of the each of the above equations as a function of the Borel mass squared, M B 2 . 
Eq. means that if the sum rule works there is a region in M B 2 where the right-hand 
side is well approximated by a linear function in M B 2 . In other words we rewrite the above 
sum rules in the form, 

g = f(M B 2 )~M B 2 -^f(M B 2 ), (30) 
dM B 

where f(M B 2 ) denotes the right-hand side of the sum rules, Eq. (129|) , and look for plateau 
in M B 2 . If we find a plateau, we conclude the value of the plateau to be the sum rule 
prediction. In this procedure the effective continuum thresholds are determined so as to 
provide the most stable plateau. 

We are now ready to discuss the numerical results of the sum rules. In the sum rules, 
we use the following parameters which determine vacuum-to-vacuum and vacuum-to-pion 
matrix elements of the quark-gluon composite operators, 

(qq) = -(225MeV) 3 , (-G 2 ) = (330MeV) 4 , 

m 2 = lGeV 2 , 5 2 = 0.2GeV 2 , e = 0.5. (31) 

We first discuss the results in the chiral limit. We searched for the effective continuum 
thresholds, uq and cu w , between 1.44 GeV and oo for which the calculated coupling constant, 
g, has the most stable plateau as a fanction of the Borel mass squared, M B 2 , over the 
range 1 GeV 2 < M B 2 < 2 GeV 2 . The optimum choices are (u , u w )= (2.4 GeV, 2.0 GeV), 
(1.9 GeV, 1.6 GeV), (2.0 GeV, 1.44 GeV) and (1.9 GeV, 1.6 GeV) for the first, second, third 
and fourth sum rules, respectively. The calculated coupling constant is plotted as a function 
of the Borel mass squared in Fig. 1. In all four sum rules stable plateaus are found with the 
value between 7 and 9.5. The OPE of the second sum rule is the same as that for the Dirac 
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structure ^5<J lJiU p fl k l/ in Ref||, while the physical parametrization of the former is different 
from that of the latter since the single pole term vanishes in the former but exists in the 
latter. However, the prediction of the former, g — 9, is numerically close to the latter since 
the coefficient of the single pole term is relatively small as is pointed out in Ref 0. The 
OPE of the first sum rule is a combination of the 275 and 275^: structures and its prediction, 
g = 9.5, is closer to the empirical value, g emp . = 13.4, than other sum rules. The third and 
the fourth sum rules do not correspond to sum rules in previous works, since the former uses 
the odd part of the vacuum-to-vacuum sum rule while the latter uses the even one. The 
calculated coupling constants of the third and the fourth sum rules are smaller than those 
of the first and the second sum rules. 

We turn to the results including 0(m n ) corrections. We first examine if the effect of the 
contninuum is regarded as correction or qualitative change in four sum rules. Figs. 2 and 
3 show g vs. Borel mass squared for ujq = ou n — 1.44 GeV and ojq = u w = 00, respectively. 
From Figs. 2 and 3 one sees that the third sum rule is very sensitive to the effective threshold 
while other three depend rather moderately on the effective threshold. Thus, we exclude the 
third sum rule from the following analysis. Then, as in the chiral limit we searched for the 
effective continuum thresholds, Uq and u n , between 1.44 GeV and 00 for which the calculated 
coupling constant, g, has the most stable plateau as a fanction of the Borel mass squared, 
Mb 2 , over the range 1 GeV 2 < M# 2 < 2 GeV 2 . The calculated coupling constants in the 
three sum rules are plotted in Fig. 4. From Figs. 1 and 4 one sees that 0(m n ) correction is 
large for the fourth sum rule which seems to indicate that the sum rule is unreliable. In the 
second sum rule the coupling constant changes almost 40% from M# 2 =l GeV 2 to 2 GeV 2 . 
Thus, the second sum rule does not satisfy the Borel stability. In contrast the first sum 
rule has a very stable plateau. Therefore, we conclude that the first sum rule is the best for 
extracting the 7rNN coupling constant. Comparing the solid line in Fig. 4 and in Fig. 1 we 
find the 0{m. K ) correction for the coupling constant in the first sum rule is about 5%. 

Finally, we disscuss the errors due to uncertainties of the parameters. We found that 
the coupling constant is most sensitive to the quark condensate, typically quoted error, 
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(qq) = —(225 ± 25MeV) 3 , changes g by ±1.7. The errors of (^-G 2 ) are expected to be 
irrelevant since the contribution from the gluon condensate is well known to be small in the 
nucleon sum rule. We next study the errors due to m 2 , and 5 2 . The parameter, m 2 ,, lies in 
the range between 0.6GeV 2 and 1.4GeV 2 according to Ref. [FJJ]. The numerical estimation, 
5 2 = 0.2 ± 0.02GeV 2 , was performed by Novikov et al. through the QCD sum rules in 
Ref. |15| . The errors due to uncertainty of m 2 , and 5 2 turn out to be less than 1% and 5%, 
respectively. We expect that the uncertainty of e is also unimportant, since from Eq. (|25"D 
its contribution is too small. We finally investigate the errors due to the uncertainties of the 
effective continuum thresholds. From the Borel stability analysis we determined the errors 
of the effective continuum thresholds as ujq = 2.5 ±0.14 GeV and = 2.0 ±0.14 GeV. These 
errors change g by ±1.3 and ±1.4, respectively. Taking account of all these uncertainties, 
we conclude 

# = 10 ±3, (32) 

which is in reasonable agreement with the emprical value. 

A comment for the choice of 5 2 is in order here. Birse and Krippa estimated the un- 
certainty of 5 2 as 0.2GeV 2 < 5 2 < 0.45GeV 2 and took 5 2 = 0.35GeV 2 in the calculation of 
g §. If we take S 2 = 0.35GeV 2 , from Eq. @ we obtain 

g « 14. (33) 

Thus the prediction of our sum rule does not become bad even if we take 5 2 much larger 
than the typical value 0.2GeV 2 . 



V. SUMMARY 

We have proposed a new approach in order to construct QCD sum rules for the 7rNN 
coupling constant, g, starting from the vacuum-to-pion correlation function of the interpo- 
lating fields of two nucleons and taking its matrix element with respect to nucleon spinors. 
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This is based on the detailed study of the physical content of the correlation function without 
referring to the effective lagrangian. In the chiral limit either the double-pole term vanishes 
or the single-pole term as well as the continuum term is not separately well defined under the 
dispersion integral in previous approaches using invariant correlation functions, while the 
double-pole term survives but the single-pole term and the continuum term give well-defined 
contributions in our new approach with projected correlation functions. In this respect the 
use of the projected correlation function seems advantageous. Another advantage of our 
approach is that unwanted pole contribution is projected out. 

Calculating the Wilson coefficients of the OPE of the correlation function up to 
O (Mb~ A ^ an d O (m^) where M B and m n are the Borel mass and the pion mass, respectively, 
we have constructed new QCD sum rules for the nNN coupling constant from the projected 
correlation function with consistently including 0(m n ) corrections. Since there are even and 
odd terms both for the vacuum-to-pion and the vacuum-to-vacuum sum rules, by dividing 
former by the latter in order to eliminate the coupling strength of the interpolating field to 
the nucleon we have derived four sum rules. 

We examined if the pole contribution dominates the continuum contribution, if 0(m n ) 
correction is reasonably small and if the result is stable in the sense of Borel stability. We 
concluded that among four sum rules the one obtained by dividing the even vacuum-to- 
pion sum rule by the even vacuum-to- vacuum sum rule is most pertinent. After roughly 
estimating errors we obtaind, g = 10 ± 3, as a result of this sum rule, which is in reasonable 
agreement with the empirical value. . We also found that the 0(m n ) correction is about 
5%. 

The present approach can be used for other interesting applications. First of all, if we 
keep the pion momentum finite in the present formalism, we will obtain the information on 
the 7riViV form factors. It is also possible to calculate not only the coupling constants of 
meson and two same baryons like rjNN and 7rSS but also the coupling constants of meson 
and two different baryons like KNU and 7rA£ without taking the flavor SU(3) limit. 
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FIG. 1. The calculated irNN coupling constant in the chiral limit with k = vs. the Borel 
mass squared, Mb 2 - The first, second, third and fourth sum rules in Eq. (29) are for solid, dashed, 
dot-dashed and dotted lines, respectively. The effective continuum thresholds, (ujq, uj t ), determined 
by the Borel stability analysis are (2.4GeV, 2.0GeV), (1.9GeV, 1.6GeV), (2.0GeV, 1.44GeV) and 
(1.9GeV, 1.6GeV) for solid, dashed, dot-dashed and dotted lines, respectively. 
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FIG. 2. The calculated ttNN coupling constant vs. the Borel mass squared, Mb 2 . The effective 
continuum thresholds are ojq = uj n = 1.44GeV for all lines. The first, second, third and fourth sum 
rules in Eq. (29) are for solid, dashed, dot-dashed and dotted lines, respectively. 
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FIG. 3. The calculated ttNN coupling constant vs. the Borel mass squared, Mb 2 . The effective 
continuum thresholds are cjo = w f = oo for all lines. The first, second, third and fourth sum rules 
in Eq. (29) are for solid, dashed, dot-dashed and dotted lines, respectively. 
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FIG. 4. The calculated ttNN coupling constant vs. the Borel mass squared, Mb 2 . The 
first, second and fourth sum rules in Eq. (29) are for solid, dashed and dotted lines, respec- 
tively. The effective continuum thresholds, (wo,w f ), determined by the Borel stability analysis are 
(2.5GeV, 2.0GeV), (1.44GeV, 1.44GeV) and (1.44GeV, 1.44GeV) for solid, dashed and dotted lines, 
respectively. 
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